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Abstract 

If a homogeneous plane light-like shell collides head-on with a ho- 
mogeneous plane electromagnetic shock wave having a step-function 



ff^ , profile then no backscattered gravitational waves are produced. We 

0|0 ' demonstrate, by explicit calculation, that if the matter is accompanied 

by a homogeneous plane electromagnetic shock wave with a step- 



function profile then backscattered gravitational waves appear after 



00 ' the collision. 
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1 Introduction 

This paper is concerned with the production of gravitational waves from 
the interaction of matter with hght in the context of Einstein-Maxwell classi- 
cal field theory. We present an explicit physical mechanism for the generation 
of gravitational waves by such a process. 

In general matter interacting with light does not produce gravitational 
waves. A simple illustration of this is the head-on collision of a homogeneous 
plane light-like shell with a homogeneous plane electromagnetic shock wave 
having a profile described by the Heaviside step-function (see p] (eqs.(3.9)- 
(3.11)) where the electromagnetic shock wave is also accompanied by a homo- 
geneous plane light-like shell with no resulting backscattered gravitational 
waves after the collision). This solution of the vacuum Einstein-Maxwell 
field equations fits the pattern of all known post-collision space-times (see 
for example [2], p]): The line-element of the space-time after the head-on 
collision of homogeneous, linearly polarized, plane light-like signals has the 
Rosen-Szekeres form [1][3] 

ds^ = -e~^ (e^ dx^ + e-^ dy^) + 2 e"^^ du dv , (1.1) 

with U, V, M functions of {u, v). The Maxwell field has only two (real-valued) 
Newman-Penrose components 0o and (j)2 which are both functions of {u,v). 
Writing (po = e~^^'^(j)o and 02 = e~'^/^(/)2, the vacuum Maxwell field equations 
read 

^ = -2^02, -Q^ = -2^u<Po, (1.2) 

with subscripts denoting partial derivatives where convenient. Einstein's field 
equations with the electromagnetic field here as source read [3]: 

(1.3) 
(1.4) 
(1.5) 
(1.6) 
(1.7) 

Consider now the head-on collision of a plane light-like shell of matter la- 
belled by a real parameter k (in the sense that ii k = then the light-like 
shell is removed) with an electromagnetic shock wave with amplitude b and 
having a Heaviside step-function profile. We take the history of the signal 
labelled by k to be found in the region u > 0, f < of the space-time with 
line-element (1.1) while the history of the signal labelled by b is found in the 



Uuv -- 


= UuU,, 


2U^u -- 


= Ul + V^-2U^M^ + Acg, 


2f/™ = 


= Ul + V:-2U,M, + A4>l, 


2K. = 


= f/„K + f/.K + 40002, 


2M„, = 


= KK-f/«f/. . 



region u < 0, f > of the space-time (the region u < 0, f < is taken to 
be flat so that the signals are non-interacting before collision). To solve the 
field equations above in the region u > 0, f > of the space-time after the 
collision requires the following conditions on the boundary of this region of 
space-time: for m > 0, f = we require 

e"^ = (1 - kuf , V = , M = 0, </)2 = 0, (1.8) 

and for f > 0, m = we require 

The solution is given by [1] 

M-ic/ _ 1-ku 

(l + 6V)3/2 ' ^'-''^ 



b 



[1 + 62^2)3/2 



(1.12) 



together with ^ = and 02 = 0. The Newman-Penrose components of the 
Weyl tensor calculated with the metric tensor given via (1.1) are 

^0 = -^(K.-t/,K + M,K) , (1.13) 

^1 = 0, (1.14) 

^2 = \{VuV,-UuU,) , (1.15) 

*3 = 0, (1.16) 

^4 = -liVuu-UuVu + MuVu) . (1.17) 

Clearly all of these components with the exception of \E'2 vanish for the so- 
lution given immediately above (because ^ = 0). This is a Petrov Type D 
Weyl tensor and since \E'o and \&4 vanish there is no backscattered gravita- 
tional radiation present. The purpose of this paper is to establish the existence 
of backscattered gravitational waves if the light-like shell here is accompanied 
by an electromagnetic shock wave. To achieve this it is sufficient to consider 
a small amplitude accompanying electromagnetic shock wave as a pertur- 
bation of the space-time with the functions U, M, 0o given by (1.10)-(1.12) 
with V = and 02 = 0. Without making the small amplitude assump- 
tion the collision problem posed here appears intractable (it is a non-trivial 



problem even with the small amplitude assumption, as will appear below) 
notwithstanding the fact that from a physical point of view it is very clear 
and simple. 

2 The Perturbed Space-Time 

If the light-like shell above is accompanied by a homogeneous plane elec- 
tromagnetic shock wave with a step-function profile and amplitude labelled 
by a then the boundary conditions (1.8) must be replaced by: ioi u > 0,v = 
we require 

e-^ = ^f^^ ,V = 0, e^ = 1 + aV , 0, = — ^ . (2.1) 

Taking the amplitude a to be small we shall neglect squares and higher powers 
of a and thus replace (2.1) by: for -u > 0, f = we require 

e"^ = (1 - kuf , V = , M = 0, 02 = a. (2.2) 

The field equations and the boundary conditions (1.9) and (2.2) give the 
following, which are useful later (in particular for determining the boundary 
conditions to be satisfied by the function /C in (2.19)): for m > 0,f = 0: 

2abu b , , 

and for f > 0,^ = 0: 

Vu = 2abv , V^ = , (j)2 = a . (2.4) 

It is convenient to write V in the form 

V^ = log([^) , (2.5) 

with A{u,v) small of order a (which we write as A = 0(a)). Thus we have 

K = 2 A„ = 0{a) , V, = 2A, = 0{a) . (2.6) 

We assume in addition that 02 = 0{a). Neglecting 0(a^ )-terms we have 
henceforth 

e-^ = (l-A;.)^ + ^-pi^-l. (2.7) 



Equations (1.3) and (1.7) give, neglecting 0(a^)-terms, 

(M + \UU = « e-"-i" = j^^^^ , (2.8) 

using the boundary conditions (1.9) and (2.2). Now the first of Maxwell's 
equations in (1.2) with 0(a^)-terms neglected yields 

- = , (2.9) 



du 
and solving this using the boundary conditions (1.9) results in 

, b 



(1 + 52^2)3/2 ' 



(2.10) 



neglecting 0(a^)-terms. Now (1.4) and (1.5) simplify, neglecting 0(a^)- 
terms, to 

2f/„„ = f/„'-2f/„M, , (2.11) 

2 U,, = Ul-2U,M, + A(g, (2.12) 

respectively. These are automatically satisfied by U,M,<pQ given by (2.7), 
(2.8) and (2.10). The remaining field equations, the second of (1.2) along with 
(1.6), will determine the two remaining unknown functions V (or equivalently 
A) and 02- 

To obtain V we begin by eliminating 02 from (1.6). We start by writing 
(1.6) in the form 

2 e-^K. = - (e-^)^ K - (e^^)^ K + 4 0o 02 • (2.13) 

Differentiating this with respect to v and multiplying by 0o yields 

00 {2 e-^K.. + 3 (e-^)^ K. + (e^^)^ K. + (e^^)^^ K} 
= -2 0^K + ^{2e-^K.+ (e-^)„K + (e-^)^K} . (2.14) 
With U and 0o given by (2.7) and (2.10) respectively we find that 

^»(^-")„,+2^2-(e-")„§=0. (2.15) 



and so the three terms in fl2.14p having \4 as a factor disappear and fl2.14p 
becomes an equation for W = V^ given by 

{{1 - kuYil + b^v^) - b\^}Wu. = -3b^v {{l-kuY-l} Wu 
+k (1 - ku) (1 + b'^v'^) W^ + 3kb\{l-ku)W . (2.16) 

To solve this equation it is useful to change the dependent variable W to a new 
variable /C(m, v) (say) in such a way that the resulting differential equation 
for /C does not have an undifferentiated /C-term. Multiplying (12.161) by v and 
writing 

"' - T^ • (^■") 

for some /C, the equation (I2.16p implies 

\{l-ku){l + 6V) - -r^^] Wu-k{l + bW) W = ICu. (2.18) 
Substituting for W from (2.17) into (2.18) results in 

h Q) Jc 

f^uv = T^ — ; — ^^v + }(^u, (2-19) 

(1 — ku)-^ 

where 

'' = " + "'"'(^-(13^) ■ <'■'"> 

The important difference between this equation and (2.16) is that in this 
equation the dependent variable /C does not appear undifferentiated (whereas 
an undifferentiated W appears in (2.16)). This makes it easier to solve (2.19) 
than (2.16). 



3 Construction of a Candidate Solution JC 

We look for the solution /C of (2.19) satisfying the following boundary 
conditions (which follow from (2.3), (2.4), (2.17) and (2.18)) : When u = 
we must have /C^ = 2abv and /Cj, = and when w = we require /C^ = 
and JCv = 2abu. We will henceforth drop the factor 2ab for the moment 
and reinstate it at the very end. 

First we rewrite (2.19) to read 

«C™-/C„ = *v{j^-i^/C„-(l-i^-L-^)/cJ . (3.1) 



This, and the boundary conditions, suggest that we look for a solution which 
is a power series in powers of 6^ of the form 

/C = Mt; + 62/C«+6"/C(2) + ... . (3.2) 

The first term here is reminiscent of a corresponding term in the Bell-Szekeres 
[B] solution of the Einstein-Maxwell vacuum field equations. The precise con- 
nection with the Bell-Szekeres solution is mentioned in section 6 below. Sub- 
stitution of (3.2) into (3.1) and equating powers of 6^ on both sides leads to a 
sequence of differential equations for K,^^\ /C*-^-*, . . .. The boundary conditions 
above result in the first derivatives of these functions having to vanish when 
ti = and when v = 0. This will fix the functions uniquely up to a constant 
in each case. We can take this constant to vanish because ultimately it can 
be removed by rescaling the coordinates x and y in (1.1). The differential 
equations for the coefficients in (3.2) read 

vK^-Ki^ = -«'+(5-^^. (3.3) 

for 2 = 2, 3, 4, . . . .These equations are straightforward to solve subject to the 
boundary conditions given above. The first few functions are given by 

^'^ - y{-"+2i(r^)' (^-^^ 

W3) 5wM 3 9 5 1 



16 [ 2k{l-ku)^ 8k{l-ku)^ 16k{l-kuYj 

(3.7) 

Substituting into (3.2) we now have 

13 5 

IC = uv b'^v^u + -b%^u b^v'^u + ... 

2 8 16 

16V / 1 \ 3 6V /I 3 2 



X] + o^ {-7;X+aX 



2 k \ 2'^J 8 k \ 2 

5 6S^ / 1 3 2 5 o , 

-7;X+oX'-T7x' +■■■ , (3i 



where the variable 



16 /c V 2^ 8^ 16 



^=0^.-^' P.9) 



has been introduced for convenience. This variable (which vanishes when 
u = 0) will appear frequently in the sequel. Remembering that we are here 
constructing a candidate exact solution of (2.19) satisfying the boundary 
conditions (the candidate will be verified to be an exact solution in the next 
section), the form of (3.8) suggests we should write 



/C 



uv 



VT+feV 



Ib'v 



2 k 

5 6S^ 



^~2^ 
1 



3bH 



,4.,5 



8 k 



^~2^+8^^ 



+ 



16 k 
1 b'^v^ 3 6^5 



X 



16 



X" 



+ 



5 b^v 



,6„,7 



2 k 8 k 16 k 
The final series in square brackets here can tentatively be written 



(3.10) 



1 



k [VTT¥^ 

Putting this into (3.10) results in 

V {1 — ku) 



1 



/C 






k^|Y^^Wv^ k 



with 



T 



2 V 2'^y8 V 2'^8'^ 

16 V 2'^8'^ 16 '^ 



From this we calculate that 



^ + 2(x + i)f 



l + 3('^)%Vx + 5g)\Vx^ 



(3.11) 



(3.12) 



(3.13) 



+7 



16 



&Vx=^ + 9(^) 6Vx^ + ...(3.14) 



We can rewrite this in the form 



^+2<'^+i)f = l;("^'- 



(3.15) 



with 

2 



g = 1 + (i) bVx + (l) bW 



+ 
This suggests that we should write 

oo 



(2n)! '' 



22«fn!)2 



{yxr = -K{yx), (3.17) 

TT 



n=0 

where y = 6^f ^ and -ft' is the complete elliptic integral of the first kind 
(see Eq.(A.l)) with argument yx- Substituting this into (3.15) we can then 
integrate with respect to x fo obtain 

T=-Ui-y,yx) + Coiy), (3.18) 

71 

where 11 is the complete elliptic integral of the third kind (see Eq.(A.3)) 
and Co is a function of integration. The boundary conditions are satisfied 
provided JF = (1 + y)^^^'^ when x = -v^ -u = 0. This is true of (3.18) 
provided the function of integration Co(i/) = 0. Now (3.12) multiplied by the 
factor 2 a 6, and with JF given by (3.18) with cq = 0, is a candidate exact 
solution of (2.19) satisfying the required boundary conditions. 



4 Verifying the Candidate Solution /C 

The candidate solution constructed in the previous section reads: 

2abv (1 — ku) Aah „, ,9 , ,, n , , , 

'= = -^7TTW ^TT" "<-""• ''"'^'- <^'^) 

It is easy to check that the first term satisfies the differential equation (2.19). 
Hence we wish to demonstrate that 

£ = t;n(-6V,6Vx) , (4.2) 

satisfies (2.19). Writing n = —iP'v'^ and m = iP'v'^x ^'^^ substituting (4.2) 
into (2.19) results in n(n, m) having to satisfy the differential equation 

d'^Ii d'^Yi dli dli 1 

2n(m-l)— ^— +2m(m-l)— ^+n^— +2(2m-l)^— + -n = 0. (4.3) 
onom om'^ on am 2 

9 



Using the formula (A. 6) for the partial derivative of 11 with respect to m this 
can be rewritten in the form 

^du , 1 . , an 1 , 

n(n — 1)— — + (— 3?72 + ?T2 + 4?T2n — 2n)— 1 — (— 3?T2 + n + 2)n 

an dm 2 

+ ^^^E-\k^O, (4,4) 

2 (m — 1) 2 

where K and E are the complete elliptic integrals of the first and second kind 
respectively (see Appendix). The formulas giving the partial derivatives of 
n with respect to m and n in terms of the complete elliptic integrals are 
given in (A. 6) and (A. 7) respectively. When these substitutions are made in 
(4.4) it follows that (4.4) is an identity. Hence the candidate solution (4.1) 
of (2.19) is indeed an exact solution of (2.19). 



5 Perturbed Fields 

In order to demonstrate explicitly that the perturbed field we are con- 
structing here contains backscattered gravitational radiation we must use the 
function /C we have obtained to calculate the leading terms (for small a) in 
^0 and ^4 given by (1.13) and (1.17) and see that they are non-vanishing. 
This we can now do because with (2.5) and (2.17), with W = V^, we have 

neglecting 0(a^)-terms, with /C given by (4.1). In integrating (2.17) with 
respect to f a function of u of integration has been put equal to zero on 
account of the fact that by (2.2) V must vanish when v = and we have 
already ensured that /C vanishes when v = 0. Since V is small of order a 
we can enter the 'background' expressions (2.7) and (2.8) for U and M in 
(1.13). This initially results in 

*»= 4"™ -2(1^ ""+!''"''"■ i^-^' 

Introducing again n = —b'^v'^ and m = iP'v'^x we obtain from (2.7) the ex- 
pression 

U, = - ,^'"-'"' ^, . (5.3) 

V [n — 1) [m — 1) 

10 



By (5.1) with /C given by (4.1) we obtain, using (A. 6) and (A. 7), 

2ab Aah f 11,, 

and 

1 3 6^^; 2a6 f (m+1) 1 

2 ™^2(l + 62^;2) ^^ - ~7rA;t;(l-A;M)(n-l)(m-l) j ^ (m - 1) J ' 

(5.5) 
with n a function of n, tti and fi', E functions of m. Thus \l/o in (5.2) has the 
non-zero value 

*° = .<.„(i-^..)(«''-i)'(m-i) '^'"-'"'"-<^"-^'" + ^'^' 

{2n — 2m — nm + 1) 3ab{n — m) 

(m-1) ^~ A;t;(l-n)5/2(,ri-l) ' ^^'^' 

which indicates that following the collision of the light-like shell, accompa- 
nied by the small amplitude electromagnetic waves labelled by a, with the 
electromagnetic waves labelled by b there exist backscattered gravitational 
waves having propagation direction in space-time d/du. In similar fashion 
we find that 

K = ---44^-.S, (5.7) 

TT (1 — KuY[m — 1) 

and 

2ahkv f ,^ (2m- 1) ^1 , , 

Tix{l-kuY{m-l)\ (m-1) J ^ ' 

The non-vanishing of this quantity indicates the existence of backscattered 
gravitational radiation after the collision having propagation direction in 
space-time d/dv. Since V is small of first order and, neglecting 0{a'^)- 
terms, U is given by (2.7), we see that \E'2 in (1.15) has its background value 
in the linear approximation given by 

kn , , 

^2 = -jz TT^^T -T^ . 5.9 

To discover what type of physical signal has as its history in space-time 
the boundary -u > 0,f = or w > 0,m = of the interaction region of the 
space-time we carry out the calculation above replacing u with u^ = u {}{u) 
and V with v+ = v 'd{v) where t? is the Heaviside step-function. We find that 



11 



in addition to the backscattered gravitational radiation found above there 
are Dirac delta-function terms in \l/o and \l/4 given by 

^^0 = - ,/^!^^ , S{v), '^, = -abv+6{u), (5.10) 

(1 — ku+) 

indicating that the boundaries of the interaction region are the histories of 
impulsive gravitational waves. We note that the Ricci tensor possesses a 
delta function term, ^Rij = —2k6{u)u^iUj, reflecting the presence of the 
light-like shell (labelled by k) having history u = 0. 

The Maxwell field after the collision has two radiative components conve- 
niently described by 0o and 02- In the approximation in which 0(a^)-terms 
are neglected we have already seen that 0o is given by (2.10). With /C and 
thus V already known in this approximation we turn to the second of the 
Maxwell equations (1.2) to find 02- Neglecting 0(a^)-terms this equation 
reads 

d4>2 b d ( }C \ 

dv 2{l + b^v^y/^ du\{l-ku)) ' ^' ^ 

Using (4.1) and (A. 3) with, for convenience, the substitution 6^f^ = y this 
simplifies to 



-/2 de 



dy 7r(l - A;-u)2(l + |/)3/2 Jo {1 - y x sin 0)^^'^ 

Noting that 



d I l-x(l + 2i/) sin2 6' \ (1 + x sin 



2 



(5.12) 



(5.13) 



dy\VTT^^l-yXsin'e) 2{l + y)V^{l-yxsm'er/^ 

we can integrate (5.12) to obtain 

02 = 2a /-/2 {l-x{l + 2y)sm'e)d0 



There is no need for a function of x of integration since evaluating this at 
V = (which corresponds to y = 0) yields the correct boundary value of 
(1)2 = a{l — ku) (see Eq.(2.2)). Rearranging the integrand in (5.14) permits 
us to write 



AaxVT+y rl'^ sin^ 9 d9 



nil- ku)"^ Jo n -U-v cln2fl^2,/l _ „ ^, cin2 



1 + X sin Oyyl — yxsin 9 

+^1 ^ ^2 /rn— n(-X,l/x) , (5.15) 

TT (1 — kuy^^l + y 

12 



with Il{n,m) given by (A. 3). Making use of (A. 7) we have 

'■"/^ sin^ 9d9 ^ (x-y)n(-x,yx) ^ K{yx) 

(l + Xsin^eyJl-yXsm^e 2x(x + l)(y + l) 2x(x+l) 



E{yx) 



2X(X + I)(t/ + 1) 



(5.16) 



and this simphfies (5.15) to read 

h = -^^m-X,yx)-K{yx)} + -^^. (5.17) 

Due to the presence of 0o and 02 after the colhsion, two systems of backscat- 
tered electromagnetic radiation exist having propagation directions d/du and 
d/dv in the interaction space-time. 



6 Discussion 

Neglecting 0(a^)-terms we have obtained explicit expressions for the 
functions U,M,V in the line-element (1.1). These are given respectively 
by (2.7), (2.8) and by (5.1) with (4.1). We also have explicit expressions for 
the Maxwell field described by the functions 0o and 02- These functions are 
found in (2.10) and (5.17) respectively. 

The approximate solution, for small a, of the Einstein-Maxwell field equa- 
tions described above has a limit k —>■ corresponding to the removal of the 
light-like shell. In this limit (4.1) gives JC —* 2abuv while \l/o -^ by (5.6), 
\l/4 — *> by (5.8) and \l/2 ^ by (5.9) so that the backscattered gravitational 
waves disappear, but the impulsive gravitational waves (5.10) remain with 
/c = in "^^o- AH of this corresponds to the important Bell-Szekeres [0] 
solution giving the exact conformally flat space-time following the collision 
of the two electromagnetic shock waves labelled by a and b. As a converse 
to the problem considered in this paper one might ask whether the intro- 
duction of a light-like shell accompanying one of these waves would produce 
backscattered gravitational radiation after collision. It is sufficient to an- 
swer this question in the affirmative by assuming small a and small k, where 
k is proportional to the energy density of the light-like shell measured by 
specified observers [2], and to specialize \&o and \l/4, given by (5.6) and (5.8) 
respectively, to this case and see that they are non-zero. We find that 



■9(2 + 76V + 46V) 7 + 76V + 36^ 



4' 



13 



^4 = -ab^kv^ , (6.2) 

neglecting 0(A;^)-ternis. In addition (5.9) becomes, for small k, 

^if2 = -b'^kv + 0{k'^) , (6.3) 

and these last three formulas confirm the limits mentioned above. 



Acknowledgment 

We thank Julien Garaud for helpful comments and Universite de Tours and 
Ambassade de France en Irlande for financial support. 



References 

[1] C. Barrabes and P. A. Hogan, Class, and Quantum Grav. 23, 5265 
(2006). 

[2] C. Barrabes and P. A. Hogan, Singular Null Hypersurfaces in General 
Relativity (World Scientific, Singapore 2003). 

[3] H. Stephani, D. Kramer, M. MacCallum, C. Hoenselaers and E. Herlt, 
Exact Solutions of Einstein's Field Equations second edition (Cambridge 
University Press, Cambridge 2003). 

[4] N. Rosen, Phys. Z. Sowjet 12, 366 (1937). 

[5] P. Szekeres, J. Math. Phys. 13, 286 (1972). 

[6] P. Bell and P. Szekeres, Gen. Rel. Grav. 5, 275 (1974). 

[7] I. S. Gradshteyn and I. M. Ryzhik, Table of Integrals, Series and Prod- 
ucts seventh edition edited by A. Jeffrey and D. Zwillinger (Academic 
Press, New York 2007), p. 859. 



14 



A Formulas for Complete Elliptic Integrals 

The complete elliptic integrals of the first, second and third kinds are given 
respectively by [7] 

/■V2 dO 

K = / - _ =K{m), (A.l) 

"'o vl — m, sm 



E = r Jl-m sin^ 9 dO = Eim) , (A.2) 

Jo 

/•'r/2 dB 

n = / , = = U{n,m) , (A.3) 

"'0 (1 — n sin 9) v 1 — rn sin 9 

where m, n are real constants. The well-known formulas for their derivatives 
are 

'"' ' 'k + ^e), (A.4) 



dm 2m \ m — 1 

'" ' 'n-H^EJ, (A.6) 



dm, 2 (n — -m) I m — 1 

9n 1 [ (n^ — m) {n — m) 



dn 2 (n — m,) (n — 1) \ n n 



n-^ -K + E\ . (A.7) 
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